Universal witness of classical correlations for two-photon states 



Karol Bartkiewicz, 1 ^ Karel Lemr, 1 Antom'n Cernoch, 1 and Jan Soubusta 1 

1 RCPTM, Joint Laboratory of Optics of Palacky University and 
Institute of Physics of Academy of Sciences of the Czech Republic, 
17. listopadu 12, 772 07 Olomouc, Czech Republic 
(Dated: February 7, 2013) 

The threshold between classical and nonclassical two-qubit states is drawn at the place when those 
states can no longer be described by classical correlations, i.e., quantum discord or entanglement 
appear. However, to check if the correlations are classical (in terms of quantum discord and entangle- 
ment) it is sufficient to witness the lack of quantum discord because its zero value implies the lack of 
entanglement. In this Rapid Communication we introduce a quantum discord witness which is able 
to detect quantum discord for a wide variety of states including symmetric states and is a universal 
classicality witness up to local unitary transformations. We also discuss how this witness can be 
practically measured in linear-optical systems using standard beam splitters and photon detectors. 
We study the efficiency of the setup assuming both ideal and real components and show that the 
efficiency of the proposed implementation is better than the full two-photon quantum tomography. 
We expect that a similar approach to the presented one can be used to directly measure quantum 
discord. The variance-like mathematical form of the discovered witness can be a manifestation of a 
nontrivial physical phenomenon possibly related to quantum phase transitions. 
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Introduction. The threshold between the classical and 
quantum world fascinates physicist since the discovery 
of quantum phenomena and realizing how different they 
are form our everyday experience. One of the prominent 
examples of quantum behavior is the non-locality lead- 
ing to violation of Bell's inequalities [TJ |5] . For two-level 
systems there is no non-locality without quantum entan- 
glement, but the opposite can be true [3]. Quantum en- 
tanglement plays a prominent role in quantum informa- 
tion processing [4] . However, the entanglement is not the 
only type of non-classical correlations. As described by 
Ollivier and Zurek [5] the non-classical correlations can 
be associated with quantum discord. Quantum discord 
(QD) is useful in many ways including quantum infor- 
mation processing or detection of quantum phase transi- 
tions, especially in the cases when the entanglement fails 
to grasp this phenomenon [5J. Moreover, it was demon- 
strated that only classical correlations can be broadcast 
locally [7] . All this features of quantum discord motivate 
the quest for developing tools for detecting and quantify- 
ing it. Nevertheless, there were only three experimental 
implementations of witnesses of classical correlations, or 
classicality witnesses (CWs), in discrete variable systems. 
Two of them were implemented in nuclear magnetic reso- 
nance systems .8, 9 and one using linear optics [lfj , how- 
ever these witnesses were not universal. At this point, we 
should stress that detecting purely classical correlations 
is a difficult problem since it involves solving optimization 
problem over not convex set of classical states. Thus, the 
problem of detecting classical correlations is harder that 
detection of entanglement. Moreover, any CW should be 
nonlinear [TT]. For those reasons the CWs [TUHTi2] are 
usually non- universal. However, there Zhang et a/.|13j 
demonstrated that finding an universal witness for clas- 



sical correlations is possible, but the established witness 
is not suitable for optical implementation. In this Rapid 
Communication we introduce CW which overcomes the 
limitations of all the previously developed witnesses. 

Let us start with introducing some basic definitions 
used throughout our Rapid Communication. A general 
two-qubit density matrix p can be expressed in the Bloch 
representation as 

1 3 
p = -(I®I+S-b t ®I+I®y-(T+'y^ j T nm a n (3a m ), (1) 

n,m— 1 

where a — [<ti, 02, 03] and matrix Ty = Tr[p((ii <8><7j)] are 
given in terms of the Pauli matrices, and Xi = Tv[p(ai £g> 
-0] {Vi = Tr[p(I ® (Tj)]) describe Bloch vector x (y) of 
the first (second) subsystem, later referred to as A and 
B. Moreover, it is always possible to transform p with 
local unitary operations p3] so that T becomes a diagonal 
matrix. 

The state p is not entangled (is separable) when it has a 
positive partial transpose, i.e., is a PPT state (see Peres- 
Horodecki criterion [mEj). The lack of entanglement 
for a two-qubit system implies, e.g., locality, in terms of 
violation of the Bell-CHSH inequality [5] (for quantitative 
study see [17]), and thus it corresponds to classical situa- 
tion where the measurement outcomes can be explained 
by a hidden-variable model. However, quantum entan- 
glement is not responsible for all the nonclassical effects. 
One of the recently celebrated manifestation of quantum- 
ness is quantum discord [SJ. The QD is responsible for 
difference in conditional quantum information calculated 
by in two ways, where one of them uses the Bayesian 
rule for calculating probabilities. Therefore, QD quan- 
tifies how much conditional quantum probabilities differ 



2 



from those calculated within classical theory. When the 
QD vanishes if the state fulfills the strong PPT condi- 
tion [18], i.e., p has to be PPT and its PPT must admit 
Cholesky decomposition (there are also other so-called 
nullity conditions - for review see [B]). Thus, if there 
is no discord, there is no entanglement. However, the 
reverse does not have to be true. 

There are several ways of quantifying QD. The one for 
which an analytic formula is known [T5j is the so-called 
geometric quantum discord (GQD) quantifying Hilbcrt- 
Schmidt distance to the closest non-discordant state. 
The expression for the GQD reads as 

D i{p) = ( Amax.t ~ ^ K,i I , (2) 

\ n=0 / 

where X n ^ (for i = A, B) stand for eigenvalues of ma- 
trix Ka = xx* + TT^ or Kg = yy< + TT\ where f 
denotes transposition. The largest A„,i is denoted as 
Amax,i- Note, that Di is asymmetric. Thus, if Da = 
the state is called classical-quantum or if Db = the 
state is quantum-classical. Naturally, there have been 
attempts of finding analytic formula for the symmet- 
ric GQD, which answers the question about the closest 
classical-classical state, however this is still an open prob- 
lem [a ng . 

If Da = Db = the state is classical-classical since 
it does not exhibit quantum correlations responsible for 
discord between conditional quantum information calcu- 
lated in the two above-mentioned ways. In the following 
sections we find an experimentally accessible criterion al- 
lowing to identify states of zero Di and describe how to 
implement it within the framework of linear-optics. 

Universal classicality witness. Here we present one of 
the most important results, i.e., the witness - Wi. It is 
constructed to fulfill Wi > 0, if the investigated states 
have nonzero GQD, and Wi — 0, if the subsystems i 
is correlated with the other one in classical way. The 
witness reads as 

Wi = - M 2 ,i > 0, (3) 

where M n ^ = Y^ m =o^mi f° r n = 1,2 are moments of 
the matrix K{ (i — A,B) from Eq. p|, where A m de- 
notes m'th eigenvalue of Ki. Note that Wi and Ki are 
asymmetric, thus Wi cannot exclusively detect classical- 
classical states. One of possible symmetric CWs is 
W s = Wa + Wb- For the witness Wi to work Ki ma- 
trix has to be positive semidefinite. 

To show that our CW fulfills the necessary condition 
for witnessing one-sided classical correlations let us start 
with expressing the GQD Di in terms of moments M n ^ 
of matrix Ki [21] . Apparently 

D L = — (A max ffil,j + A m ax02,i) (4) 
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FIG. 1. (Color online) Setup for linear-optical measurement 
of U m n (Vmn) observable. The setup consists of 50:50 beam 
splitters (BS) and conventional detectors (D). The possible 
measurement outcomes for U mn (Vmn) are -4 if D2 and D3 
click simultaneously as we witness singlet projection or 1 (2) 
if coincidence is registered Di and D4 (this measures /). The 
device works with probability 0.5. 

and Wi = M\i — M2.1 = 4A maXi i_Di — A max ®i^A max ©2,i, 
where © stands for sum modulo 3. It is apparent that 
for A m ,i > the witness Wi fulfills the necessary condi- 
tion. This is because for positive Ai the vanishing of GQD 
given in Eq. Q implies A maxeM A maxffi 2,i = 0. Thus, the 
nullity of witness Wi is necessary for Di = 0. 

To demonstrate that Wj is a sufficient witness of clas- 
sical correlations it we have to show that 

Amax©l,iA max ®2,i = A max ,i ( A max ® + A ma x©2,i) (5) 

holds only if both sides are zero. Since A max .i > A n ,i 
for n — 0,1,2 the equality is possible only if one of 
A n ,j = 0, say A max ®2,i = 0, but this implies that 
A maX ffii,iA ma xe2,i = and the only possibility for the 
equality to be satisfied is to have \ max ®i,i = 0. Since 
A max ©i,i = and A maxffi 2,* = we must have D. L = 0, 
because Di it is directly proportional to sum of those 
eigenvalues. Thus, our witness is a sufficient witness of 
quantum discord. 

The witness Wa (Wb) is applicable as long as Ka 
(Kb) is positive semidefinite. Since, xx* or yip is pos- 
itive semidefinite, the requirement for Ki reduces to re- 
quiring TT^ to be positive semidefinite. This obviously 
happens, e.g., for symmetric states pab = Pba which in- 
clude Bell-diagonal states and other states like symmet- 
rically amplitude-damped or phase-damped Bell states 
[T71 12"2"] . However, we also get Ki > if T is diagonal 
which can always be assured by local unitary transfor- 
mations on systems A and B [14]. Hence, our classicality 
witness can be considered universal up to local unitary 
operations. Moreover, since the symmetric geometric dis- 
cord vanishes D s = only if Da = Db = [20]) for 
checking if the state exhibits classical-classical correla- 
tions W s = Wa + Wb should be used. 

Linear- optical implementation. Up to this point we 
did not mention how to measure the CW. The proposed 
measurement procedure is similar to the proposal for di- 
rect measurement of quantum discord introduced by Jin 
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FIG. 2. (Color online) Setup for linear-optical measurement 
of moment (a) Mi and (b) M2. The setup works with prob- 
ability 1/4 (1/16) for measuring Mi (M 2 ). Note that the 
measurement of Mi and M2, and thus W, can be conducted 
in parallel if one has access to six copies of the input state p. 

et al.m [5T] and it resembles the so-called entanglement 
swapping [23 . 

In our case we need only first two moments of matrix 
Ki, and we need a fewer copies of p. The first moment is 
given as Mi a — Tr(xx" + TTt), where wc fixed without 
loss of generality i = A, and requires two copies of p to 
be measured. This is because M iy A can be also expressed 
[H] as 

M M = Tr[(U AlA2 <g> V Bi b 2 )(pa iBi <8> Pa 2 b 2 )}, (6) 
where 



FIG. 3. (Color online) Optical setup for measuring M%. When 
the delay lines are removed, the setup measures Mi as in 
Fig. [2^,. At time t = and t = r two photon pairs p enter 
the setup. When t — 2r, 3r no new photons appear and 
only coincidences in the left part of the setup are witnessed. 
This part consists of U measurement, where one photon is 
always delayed by r and the second one randomly by 2r or 
0. The second delay works with probability of 0.5. There is 
no delay in the right part of the setup measuring V. In the 
successful cases pa x b x and pa 2 b 2 enter the setup at t = 
and Vb x b 2 is measured. At the same time photons in the left 
part are delayed, the first from the left by r and the second 
by 2r (with probability 1/4), so at t — there is no photon 
detected in U. Photon pairs Pa 3 b 3 and pa 4 b 4 enter the setup 
at t = t and Vb 3 b 4 is measured without delay. In the left 
part, the first photon is delayed by r, but the second one is 
not affected (with probability 1/4). At the same time Ua 1 a 4 
is measured. Since there are no photons added at t = 2r 
only Ua 2 a 3 is measured (for a photon provided at t = and 
delayed by 2r and a photon provided at t = r delayed by r). 
The success rate of the procedure is 1/256. The 1/16 drop 
in comparison to Fig. [2Jd is caused by applying probabilistic 
delay. The whole sequence is repeated until a good estimate 
of M2 is obtained. No photons enter the setup at t — 3r to 
not affect the next iteration. 



V A m A n = ~ 4P A A +lA m A nl VB m B n = U Bm B n +lB m B n - 

(7) 

Both U and V are given in terms of singlet projections 
Pi min = l*->i m in<*"l and |*-) imin = (|0} lm |lk - 
|l)iJO)iJ/V2 with i = A,B and to, rt = 1,2. However, 
for measuring M^^a we need four copies p, since 

M 2 ,A = ^[{U Al A 4 ® Ua 2 A 3 ® V Bl B 2 ® Vb 3 B 4 ) 

X (fiAiBt <8> P^ 2 S 2 8) PA3B3 ® PAi-bJ]- ( 8 ) 

The projections P - and J can be implemented using 
linear optics (see Fig. [TJ). The operators V A m A n and 
Vb„ b„ are local two-qubit operators acting on the respec- 
tive subsystems. Having four copies of p each moment 
can be estimated in a single coincidence measurement. 
However, if six copies are available, the CW from Eq. ^ 
could be evaluated in single measurement. Multiphoton 
experiments are demanding but feasible (for review see 
Ref. [H]). Nevertheless, with a few modifications the 
measurement can be implemented as presented below. 

By choosing orthogonal photon polarization states as 
z-basis for subsystems constituting p we can implement 
the above-introduced witness from Eq. ^ in a linear- 
optical system. Our proposal for the experimental setup 



consists of standard optical elements, i.e., 50:50 beam 
splitters (BS) and photon detectors. The implementation 
requires up to four elements from Fig. [T] This means 
using simultaneously up to four photon pairs in state p 
as shown in Fig. [2] 

The U and V can be measured using the device shown 
in Fig. [TJ where with probability 1 /4 both photons pass 
unchanged through the pair of the 50:50 BS. Next, if both 
photons are detected by the D\ and D4, we measure the 
part of the operator U or V proportional to /. Alterna- 
tively, with probability 1/4, both photons reflect on the 
BSs and overlap on the central 50:50 BS. After overlap- 
ping on the BS the state p is projected onto a singlet 
state by witnessing coincidence in D2 and D3. A setup 
that can be used for measuring M^j (M2,i) consists of 
one (two) V mn and one (two) U mn element (s) arranged 
as shown in Fig. 2a(b). However, obtaining valid co- 
incidences is probabilistic and thus the success rate of 
measuring Mi is 1/4 and for M2 it reaches 1/16. In 
spite of loosing many photons measuring Wi is expected 
to be much faster than full 2-qubit quantum tomography 
(QT) which requires rotating a number of polarization 
plates providing 16 measurement configurations, where 
each rotation takes a few seconds. During this time the 
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measurement setup from Fig.[2]would accumulate enough 
data to estimate Wi- 

The main difficulty in implementing the setup out- 
lined in Fig. [2] is working with large number of photons. 
The photon pairs should exhibit quantum correlations, 
thus should be produced, e.g., in spontaneous parametric 
down-conversion of light of a certain degree of depolariza- 
tion. Since the pairs are produced by a random process, 
the coincidence rate in our experimental proposal can be 
very low. There are of course further challenges. De- 
pending on the approach of producing the input states, 
one may face the problem of indistinguishably of the pho- 
ton pairs, which causes the results to differ from the ex- 
pected one, but methods of circumventing this effect can 
be found, e.g., in [24]. Moreover, the theoretical predic- 
tion is reached if all the pairs perfectly overlap on the 
corresponding BSs. This condition is easy to be sat- 
isfied for one pair of photons, however the amount of 
work increases with the number of photons. To achieve 
this all the optical elements would require active stabi- 
lization. Nevertheless, we can reduce complexity of the 
optical setup by using temporal separation between pho- 
ton pairs p as shown in Fig. [3l The success rate of the 
measurement of Mi from Fig7|3] is 1/256 instead of 1/16 
as in the setup from Fig. since it requires applying 
probabilistic delay for measuring U. When all the de- 
lay lines in Fig. [3] are removed, the setup corresponds 
to the one shown in Fig. [2^, and can be used for direct 
measurement of M\. By taking account for the pauses 
in providing two pairs of photons at t = 2r, 3t (we only 
allow every second two pulses from the pulse train to en- 
ter the setup by using a pulse picker) the efficiency of the 
setup drops by a factor of 1/2. Thus, for perfect detec- 
tors, the useful coincidence count rate would be limited 
to 1/512 of the brightness of the source. This effectively 
changes the repetition rate of the source from 80 MHz 
to 160 kHz. Moreover, this rate is further reduced by a 
factor of 1/2 because of the 50 ns dead time of the de- 
tectors. Finally, we arrive at 80 kHz effective repetition 
rate of the source. In two-crystal type-I geometry (so- 
called Kwiat source [IS]), the probability of generating 
an entangled photon pair from a single pump pulse is 
about 1/10. Thus, creating two pairs independently oc- 
curs with probability 1/100. Assuming perfect detectors, 
we can expect the photons generated with a pulsed laser 
at repetition rate of 80 kHz to cause the 8 photons (two 
pairs generated at t = and other two at t = r) to ar- 
rive at the designated places 800 times per second. In 
reality this number would be much lower due to finite ef- 
ficiency of the detectors. For realistic detector efficiency 
of 0.75 (Perkin-Elmer single photon counting modules 
operating at wavelength of 700 nm), the number of coin- 
cidences per second would reach about 80. This number 
would be further reduced by a few percent by unavoid- 
able imperfections of the setup. Nevertheless, the final 
expected number of detection events for Mi is of order 



of tens coincidences per second which would allow us to 
measure Mi in less than a minute. Measuring M\ would 
be much faster because of the simpler structure of coinci- 
dences, no pulse picking, and lack of probabilistic delay. 
The original setup for measuring Mi outlined in Fig. [2]d 
would also not suffer from such a low coincidence rate as 
the one shown in Fig. [3} However, the later alternative is 
more feasible and can be successfully implemented, e.g., 
in our laboratory. 

Conclusions. We discovered a novel CW (witness of 
GQD) which is universal up to local choice of coordi- 
nate system. Our witness will work as long as the Ki 
matrix is positive semidefiite, which includes, e.g., sym- 
metric states like Bell-diagonal or symmetrically-damped 
states [T71I22]. Therefore, our CW is more powerful than 
the ones reported recently by others [8Tll2|. In this re- 
spect can our CW compared with the one introduced by 
Zhang et aim [13 which also requires four copies of the 
input state, but is not appropriate for optical implemen- 
tation. Our approach is less complex than full QT, since 
it only requires 2 measurements, if 4 copies of p are avail- 
able [either spatially (see Fig. [2]) or time separated (see 
Fig.(3|], instead of 16 for the QT [26H28]. Furthermore, 
the proposed measurement setup is expected to be faster 
than the direct QT, since it would provide the outcome 
within a minute. Let us note that if the one had access 
to 6 copies of p, the witness W could be estimated in 
a single measurement in contrast to the simpler experi- 
ment described in |10j . where the CW was not universal 
and its estimation involved measuring three quantities. 
Furthermore, our approach can be extended to directly 
measure GQD as discussed in Ref. [21]. We expect the 
simple mathematical expression for the CW to have a 
deeper physical meaning since it strongly resembles the 
expression for variance where role of the stochastic vari- 
able is played by eigenvalues of K t and might be related 
to quantum phase transitions. 
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